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Abstract 

We study the depinning phase transition of a directed polymer in a (i-dimensional space by 
a periodic potential localized on a straight line. We give exact formulas in all dimensions for the 
critical pinning we need to localize the polymer. We show that a bounded state can still arise 
even if, in average, the potential layer is not attractive and for diverging values of the potential 
on the repulsive sites. The phase transition is of second order. 

PACS numbers: 68.45. Gd; 68.35.Rh 
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The statistical mechanics of long linear chains (directed polymers) in ordered and disordered 
media has been object of intense study in the past years. One relevant problem in that large 
context is the study of the depinning transition of a directed line by a single extended defect 
embedded in a {d + l)-dimensional lattice. An attractive potential pins the interface (here 
a line) on itself suppressing wandering, but thermic fluctuations increase the configurational 
entropy and a phase transition takes place at a given critical temperature T c The depinning 
transition of polymers by a single defect has been object of intense work and exact results are 
now well established in some simple cases Moreover the simultaneous effect of both point 

and extended defects has lead recently to some new results in the context of RG approach in 
continuum models 

In this article we will deal with the following problem: the depinning phase transition of a 
single polymer in a (d + l)-dimensional hypercubic lattice by a periodic potential localized on 
a line, i.e. a potential which is alternatively attractive and repulsive. We point out that such 
potential layer can be used, for instance, to mimic the effect of two alternating kinds of pinning 
centers with different strengths. Moreover it is tightly linked with a simplified version of the 
KPZ equation for interface growth (see discussion below). Some well-known arguments show 
that with a n-dimensional (oriented) defect the polymer is localized for d — n + 1 < 2 by an 
arbitrarily weak attractive force, while for d — n + 1 > 2 a finite strength is necessary to do the 
work |j. Analytical results are in general not available for such high-dimensional systems. 

Our main results in this work are the following: (i) we solve the phase transition problem 
in all dimensions and we give an exact formula for the critical pinning strength necessary to 
localize the polymer at the origin, (ii) We surprisingly show that a bound state can always arise 
even though, in average, the potential layer is repulsive. This effect has been recently studied 
in the one-dimensional case ||. We also prove that in all finite dimensions a finite strength on 
the attractive sites is enough to pin the polymer even in the limit of an infinite potential on 
the repulsive sites, (iii) The approach introduced is also interesting by its own: we use a dual 
space representation of the transfer matrix which enables us to simply find the critical state 
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of the system and the partition (wave) function. This approach leads to exact results for the 
error-catastrophe problem in biological evolution ]?]]. 

The energy of a line of length L with extremes at = and = x, wandering in a 
(d + l)-dimensional space SI = Z d ® N and directed along a "time" axis, is given by Q: 

H({h^}) = jj2\^ k) -^ k - 1) \-J2 U ^,5 > (!) 

k=l k=l 

where b_( fc ) is a vector identifying the position of the line in Q at each "time" k. The potential 
is localized at the origin and it is alternatively attractive and repulsive, i.e. XJ\- = u > if k 
is even and = —v < if k is odd. The directed line has no overhangs and RSOS condition 
is imposed. A canonical partition function is introduced (the sum is over all possible allowed 
realizations of the interface "height" h^) : 

z L (x) = $> P (-w({h( fc )})/r). (2) 

{h} 

In the usual approach one defines a symmetrical transfer matrix T zz / = T z - Z > from Z as 
Ji = — TlnJ2z T(«)S Z ; (here S z stands for a shift operator jlj]). At finite temperatures T > 
the fluctuations of the interface increase the configurational entropy while large humps are 
unlikely since they give a higher internal energy. The final state of the polymer is the result 
of that competition and it is associated to the free energy density (per unit length) /. In the 
thermodynamic limit (L — ► oo) / is dominated by the largest eigenvalue of T zz '. 

For our system we see that in one step the partition function -Zx(x) obeys the following 
recursion relation: 



Zl+i(x) = [1 + (a L+ i - 1)<5 X>5 J (3) 

L 

x Z L (x) + t]T(z L ( X + e«) + iJ L (x-e«) 



L i=l 

where the unitary vectors eW = (0, 0, • • • , i, • • • , 0) have a "1" bit as i-th. element (so to 
satisfy RSOS conditions). In the above we have defined the parameters t = exp(— J/T) £ (0, 1] 
cll = exp(n/T) (for L even) and ai = exp(— v/T) (for L odd). We now introduce a dual space 
representation of our equation in order to simplify the calculation. For the present problem we 
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use a standard Fourier transform, but sometimes one needs different representations j7j. As the 
partition function is expected to be symmetric in the arguments we introduce a cosine 

transform 

Z L (x) = ( l & d k TTcoB(7rfciSi)Z L (k), (4) 

and its inverse. In the Fourier space eq. (3) takes a simple form; a proper definition of the 
transfer matrix should anyway take into account the periodicity of the problem. If we introduce 
the normalized quantities G L {k) = Z L (k)/(l+2dt) L and£(k) = (1+2* J2i=i cos(vrA; i ))/(l+2dt), 
after two consecutive steps our equation reads: 

G 2L+2 (k) = e 2 (k)G 2L (k)+A /Vg^q^q) 

J o 

+ B£(k) /Vg^G^q) (5) 
Jo 

A R rl 

+ msi Af(q)G 2i(q ). 

with A = a 2 L — 1 an d B = a 2 L+i — !• 

In order to find the maximum eigenvalue e we should consider the spectral equation obtained 
from (5) by identifying Gx_|_ 2 (k) = e 2 Gi(k). We recall that the only significant contribution 
in the thermodynamic limit L — > oo is the maximum eigenvalue associated with (5). In fact 
one can show that / ~ — loge is always different from in the localized region (e > 1) while it 
vanishes in the unbounded state for e — > 1. The depinning phase transition is defined at e = 1 
I- 

The search for a general solution of the eigenvalue equation is a very hard task, nevertheless 
one can find the criticality condition as a function of the free parameters {d, t, u, v} of the theory 
with no enormous effort. The idea is to introduce two auxiliary constants i^i j2 by integrating 
all the terms containing the unknown function Gx(k): 

K n = f 1 d d k G L (k)C(k) (n = 1, 2). (6) 
J o 

Therefore we get a homogeneous system of two algebraic equations for K n which must be 
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satisfied by any general set of parameters {d, t,u,v}: 



AB 

Al 2 {e) + T ^-li{e)-l)K l +Bl l (e)K 2 = 0, 



(M 3 (e) + Y ^l 2 (Ef)K 1 + (Bl 2 (e)-l)K 2 = 0, (7) 



l n {e)= d d k 2 r( k) (n = l,2,3). (8) 



with 

T - [ A d h 

e 2 -£ 2 (k) 

The homogeneous system (7) admits non trivial solutions iff the determinant of its coefficients 
vanishes. This is therefore the condition we must require in order to get the spectrum of the 
transfer matrix. Performing the limit £->lwe arrive at the condition which must be satisfied 
by any set {d, t, u, v} at the critical point [||]. After some calculations we then get the criticality 
condition: 

1 - (A + B)l> 2 + AB [z* -2< (l' 3 + — ^)] = 0, (9) 

with X' n = T n {e = 1), (n = 1,2,3). In the following we will drop the "prime" from the 
formulas, anyway recalling that all quantities are calculated at e = 1. Borrowing from the 
thermodynamic language, equation (9) can be thought as the equation of state at criticality: 
the "thermodynamic variables" are now those in the set {d, t, u, v}. 

Despite of the complexity of the high dimensional integrals 2 n involved in the above formula 
one can finally express them, after some analytical work, in the following form: 

l 1 = ^(f-g), J 2 = |(/ + 5 )-l, T 3= J 1 __L ) (10) 

where we have defined a new constant a = (1 + 2dt)/2t and the two integrals: 

with J2' = J2i=i cos(irki). By means of the above definitions our equation of state reads 



Bl 2 -1 
B{ll-ll)-l 2 



a = r ■ (12) 



The integrals / and g are well known in the theory of random walks (RW) Q; the former, in 
particular, has a well defined physical meaning: it gives the mean time spent on the origin for a 
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random walker in a d-dimensional hypercubic lattice (times 1/d). In other words 1 — d/f is the 
probability of return of a random walker to his starting point. We use an integral representation 
to write them in a simpler form: 

/■oo roc 

f= / due- du I (u) d , g= / du e-( d+l ^ u I (u) d . (13) 
Jo Jo 

Here Io(u) is the usual modified Bessel function of integer order. The mathematical properties 

of / and g are central to our solution and then we will summarize them in more detail. All 

below results hold in the ranges: d G [0, oo) and t G (0, 1] ||]. 

Both / and g are positive strictly convex decreasing functions of d, converging to + for 

d — ► oo. Moreover we have that / > g \/t G (0, 1]. It is interesting to look at their behavior in 

some extreme situations. If we perform an asymptotic development of g for t close to we get 

the result g = t — dt 2 + 0(i 3 ). Moreover we find that 

f = ?(l + l + ^ + ^S + o(l)) (14) 



d V 2d Ad 2 2d? \d A , 

1 d 3d 2 n (l_\ 

9 ~ d + 1/t + 2(d + l/tf + 4(d + l/tf Uv 

from the asymptotic developments at large d. Perhaps the most important property of the two 
above integrals is that / is a divergent integral for d = 1, 2 while it is finite for d > 2 (e.g. see 
II). The second one, g, is finite, on the other hand, Vd. In the RW theory the divergence of / 
for d < 3 leads to the well-known result that the total probability or return of a random walker 
to his starting point is 1 only for dimensions less than 3. Our solution shows directly how this 
pure topological effect plays a central role in the context of the depinning transition for directed 
polymers. 

For d = 1, g can be explicitly calculated and by taking the dominant contribution of (12) at 
diverging /, the criticality condition becomes 

S=l°g(l 7 ~ - \ 1 , (16) 



which confirms the result obtained by Nechaev and Zhang in the same context ||. Let us now 
turn back to eq. (12) for the general case. At t = 0, or equivalently J — > oo, the polymer 
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is a rigid straight line and then it is in the pinned (resp. unpinned) phase depending on the 
sign of the difference u — v. This result is indeed contained in our solution: by asymptotically 
expanding g for small t (see above) we finally find 

B 4 + 4B-B 2 (df - 1) n . o, 

We see that, as expected, at vanishing hopping constant t there is no more dependence on the 
dimensionality; and solving for the pinning strength one gets the result that A(B + 1) = —B 
or u = v. The "phase diagram" on the plane u-v is then represented in this case by a single 
straight line bisecting the whole space (see Fig.(l)). In the above semispace (u > v) the polymer 
is in a bound state, while in the lower one (u < v) it is completely delocalized. 

What does happen if the hopping constant t is different from 0? In, by now, standard 
notation we define u c as the value of the force on the attractive sites which satisfies (12) for a 
given set {d,t,v}. For B = 0, or equivalently v = 0, the potential layer is made of alternating 
attractive and neutral sites and the criticality condition reads u c = Tlog(l + I/Z2). This is an 
exact formula valid in all dimensions. Since / (and then I2) diverges for d < 3, an arbitrarily 
small u = 5 > is enough to localize the polymer, according to well-know general results, while 
for d > 3 we need a finite value 
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u IT OJ -H2t-2t 2 l-2t + t 3 -3t 4 l „ 
6 ' = M+ * + 2^ d + ° 



(18) 



t 6 d 2 _ 

to do the job (this phenomenon has also a QM counterpart [0]). Moreover that critical attrac- 
tive force diverges logaritmically at large d. The same result we get in the more familiar case 
v = — u representing an extended linear attractive defect. In this case one finds 

1 



1 + 



(19) 



X1+J2. 

which gives, at d = 1, a known result, i.e. u c = 0. Fig. (2) shows the shape of u c /T as a function 
of the dimension d in that case. 

An interesting aspect of our solution is that a bound state can take place at all finite 
dimensions d even if the potential layer is, in average, repulsive (v > u). This is evident from 
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Fig.(l) in which we have drawn the critical curves separating bounded and unbounded regions 
for different values of d as a function of the "reduced parameter" v/T. Above (below) the 
straight line u = v the potential is, in average, attractive (repulsive). Then we see that for 
every fixed v and VcZ (finite), one can find a finite value of u c < v which localize the polymer 
giving a bound state for the partition function. The more astonishing point is that one can 
simply prove from (12) that the critical curves asymptotically converge toward a finite u c also 
in the extreme limit v — > +00: 



(00) 
Uc . 



1+ 1+12 



(20) 



(J|-T 1 2 )+T 2 . 

Again, for d = 1, we recover the known result §] exp(n^ oo) / T ) = v / T+~27/(v / r+~2T - 1). The 
intuitive explication of this apparent paradox is that the polymer wanders in the space avoiding 
repulsive sites and passing through the potential layer on the attractive ones (preferentially). 
Some critical curves are showed in Fig.(l) for both d < 2 and d > 2 as functions of the reduced 
parameter v/T. 

We recall that, as obvious, instead of choosing u and v as "free parameters" in the above 
considerations, we could, in principle, directly look at the behavior of the system as a function 
of the temperature T (from eq.(12)) with u and v fixed. 

In presence of a potential layer the broken translational symmetry of f2 can be associated to 
the presence of Goldstone modes whose mass fi is finite for T < T c (i.e. in the pinned phase). 
The bounded state can be expressed as 2(x) = exp(— /x|x|) and the maximum eigenvalue as 
e = 1 + /i 2 . By using the explicit form of the partition (wave) function in the localized phase Q 
we find that near the transition \ f — f c \ vanishes quadratically with (u — u c ), i.e. the transition 
is of second order. Then the transversal correlation length £j_ diverges as (u — u c ) _1 at the 
critical point. 

Conclusion: we have studied the problem of a polymer in a d-dimensional space in presence 
of a linear extended defect with periodically arranged pinning sites and we have found the exact 
condition for occurring of the depinning transition. Exact formulas for the critical pinning are 
found in all dimensions. Perhaps the most interesting aspect of our system is that a bounded 
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state of the partition function can arise for all d even if the potential on the extended defect 
is in average repulsive. We could also get a deeper insight into the physical problem from the 
calculation of the of the partition wave function. This seems to be a very interesting aspect 
which can also be dealt with our approach || . Recently has been pointed out that the physics of 
a polymer subjected to a linear extended defect can be seen, in the context of dynamical growing 



of interfaces, as a simple version of the KPZ equation [11] with a "noise" term proportional to 
8 (x) | |12| |. We will present our results on these problems in a separate paper. 
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Figure captions 

Fig. 1 

Critical curves u c /T vs. v/T for different values of t and d calculated by numerically integrating 
/ and g (see text). Above (below) them the system is in a bounded (unbounded) state for the 
partition function. At vanishing hopping constant (t=0) the depinning line is the diagonal of 
the phase space u-v. For t > and d > 2 we need a finite value of the potential u to pin the 
polymer. 

Fig. 2 

The critical pinning (divided by T) necessary to localize the polymer as a function of d in 
the case of a uniform attractive potential (i.e. v = —u). The divergence of u c /T with d is 
logarithmic (see text). 
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